B&P File No : 8389-018 
BERESKIN ft PARR 



UHITEP STATES 



Title; Development of a Compact Raman Spectrometer for Detecting 
Product Interfaces in a Flow Path 

Inventors,: H. A. Gamble, J. C. Robbine, G. I. Mackay and H. I. Schiff 




Title: Development of a Compact Raman Spectrometer for Detecting 
Product Interfaces In a Flow Path 

5 

D OF THE INVENTION 

The present invention relates to a method and apparatus for 
detecting transitions between different gas or liquid products in a tlow path 
and, more particularly, it relates to an appaialus and method utilizing Raman 
10 spectroscopy for detecting transitions between petroleum products. 

BACKGROUND OF THE INVENTION 

Before the advent of lasers, the use of Raman spectroscopy as a 
routine analytical technique was limited due to the small number of available 
15 sources of intense, monochromatic radiation Since the 19B0's, lasers have 
become tne excitation source of choice fui Raman spectroscopy, as they 
provide rnuuh yiealei inturibillub and naiiuwei liuw widths than the mercury 

arc lamps commonly employed previously. Furthermore, much weaker 
Raman signals became observable, resolution improved which lead to laser 

20 Raman spectroscopy becoming an important benchtop tool for identifying 
molecular species via characteristic or "fingerprint" wibrationa 1 features. In the 
past, due to the special requirements of the lasers, (high voltages, cooling 
water, specialised personnel and space requirements), tnese systems tended 
to be large and expensive, and needed to be used in dedicated facilities. 

25 Since then, diode lasers have become much simpler to operate. Diode 

lasers are small and inexpensive, can run off very low vo!tages (15 V or less), 
geneiale less heat, and nave high conversion efficiencies compared to 
traditional laser Raman sources, in spite of these advantages, certain inherent 
properties of diode lasers have made them less appealing for use in Raman 

30 spectroscopy. These include lower intensities, a less monochromatic output 
( mode hopping"), greater beam divergence, and excitation wavelengths 
restricted to the near to mld-lR regions. With recent advances in diode lase* 
technology, many of these difficulties have been overcome, and Raman 
spectrometers with diode lasers as the excitation source have begun to 

35 appear 
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One big advantage of using diode lasers as a source for Ramar 
spectroscopy was to move the technique out of the Ian and into the realm Of 
fieio measurements and process monitoring, Relatively high molecular weight 

5 organics, such as petroleum products, various plastics, and types of edible 
oils have oeen differentiated on the basis of their Raman spectra generated 
using a portable diode laser based Raman instrument. These compounds are 
good candidates for Raman based analysis due to the presence of strong 
Raman features with shifts in the region of 700-1700 cm' 1 . Raman analysis of 

10 commercial grade gasoline in particular can benefit from an excitation source 
in the near to mid-IR, to avoid interference from background fluorescence 
which can be excited at lower wavelengths. These fluorescence signals 
depend on the excitation wavelength, and can be so strong that they 
completely ooscure the Raman features which would otherwise appear. 

15 A number of studies have used Raman spectroscopy to examine fuels 

0' mock fuel mixtures. These include a quantitative analysis of xylene isomers 
in mock petroleum fuels using a diode laser Raman spectrometer with an 
epilation wavelength of 800nm. A partial least squares regression analysis 
routine was used to correlate the individual xylene isomer concentrations to 

?o the Raman signal, without the use of an internal standard. For samples 
containing between 1.5 and 15% xylenes, the concentrations were 
determined to within ±0.1% for the rt*tm- and pa/-a- isomers, and la within 
±0.16% for ortto-xylene. Other studies include a comparison of near-lR and 
Raman spectroscopies for the determination of the chemical and physical 

25 properties of naptha an analysis of aviation turbine fuel composition, and a 
system designed lo correlate the Raman spectra of gasolines with their 
octane ratings. The first three studies were all laboratory based, while the 
fourth describes a partial least squares regression analysis routine which was 
applied to spectra recuided on a commercial FT-Raman spectrometer with a 

30 Nd: YAG source. A large "training eef of spectra taken from fuels with known 
nctane ratings was used to build a mode! to predict the octane rating of fuels 
not included in the training set. The accuracy of the determined octane ratings 

depended on the accuracy uF lite- Liaiiiiiij se.l usftd to create the model In 
ycucial, o given fuel octane rating can correspond to any number of different 
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chemical 'recipes', i.e the octane rating does not uniquely deHne the exact 
chemical composition of the fuel. Casoline derived from a common source, or 
procaeeod by a particular refiner/, may show a particular pattern which the 
Training set can ■•learn" to recognize. However, a fuel derived from a different 

5 source may not match the defined pattern, necessitating the acquisition and 
use of a new training set. 

The problem of measuring the octane rating of a given fuel, other than 
by the empirical "engine knock" tests used to define the quantity, is clearly not 
straightforward l he training sat approach does have Its uses. However, it 

10 relies on berg able to establish a reliable base (the training set), and can get 
somewhat cumbersome when applied to a large, widely varying set of 
samples. The probtem of distinguishing between various grades of gasoline, 
without attempting to specify an octane rating, is significantly less demanding 
It is also still extremely important. When the finished products leave the 

15 refinery, they travel through pipelines to distribution stations, where they are 
directed into holding tanks before being transported by truck to local filling 
stations. Any or ail of the different grades of gasoiine or distillates may pass 
sequentially thiough a given section of pipe. It is thus important to be able to 
determine exactly when one product lot ends, and the next begins. 

2n Vanous techniques have been used to identify the exact product 

interface in such a setting. As a basic requirement one needs a measurable 
physical or chemical property, which differs not only from product to product, 
but from grade to grade. Ideally, the measurements should be fast, non- 
destructive, able to be made in-srtu, give a clearly visible (large) change when 

25 an interface occurs, and produce an output, which does not require a high 
degree of technical skill to interpret. 

One technique is to measure the density of the products ae they flow 
past a particular point near the outlet into the holding tanks The density 
measured al Ihis point in actually based on a sonometer, which measures the 

30 speed of sound as it passes through the sample. This technique produces 
large changes at gasoline to distillate interfaces, but can have trouble 
distinguishing between grades of gasoline, and between diesel vs. low 
eulphur diese! oils Viscosity or colour changes have also been used Another 



device, advertised commercially, detects inrerfaces by measuring the 
electrical resistivity of Hie product as It flows past a point. 

The abundance of techniques available serves to demonstrate the 
importance ot detecting product interfaces in a gasoline pipeline. Accordingly, 
* tnprp is » need for a method, which reliably and accurately detects product 
interlaces in a pipeline in real-time. Particularly, there is need for a method 
and apparatus which can accurately and reliably detect product interfaces 
between a wide range of petroleum products, such gg gasoline products as 
well as cthe* distillates (Diesel, jet fuel, etc ). 

10 

SUMMARY OF THE INVENTION 

In accordance with the present invention a compact, portable Raman 
spectrometer system suitable for in situ use in hostile environments has been 
proposed. The source is a high power (500 mW), broad band red diode laser 

1b which has been mode locked using an externa! cavity. This produces a 
monochromatic excitation beam at a wavelength of approximately 670 nm. 
The spectrometer consists of an entrance siit, a combined diffraction 
grating/focussing element, and an exit slit. The resolution of the Raman 
p.nenrra nhtained is excellent (7.2 cm' 1 FWHM at a Raman shift of 1000 cm" 1 ). 

20 The Raman signal, which exits the spectrometer exit slit is detected by a 
highly sensitive photomultiplier tube (PMT), and sent to a computer device 
(PC-1C4) for data acquisition and analysis. 

The proposed invention described herein detects liquid or gas products 
in a flow path. Specifically, this invention detects changes in the composition 

25 of various petroleum products flowing through a gasoline pipeline, by means 
of exposing samples of various petroleum products to the Raman 
spectrometer system. The petroleum products of interest consisted of four 
distillates (diesel oil, low sulphur diesel, jet fuel and furnace ml), and three 
different grades uf gasoline (reyular unleaded, premium 91 and premium 

30 US92). Both the hardware and software is tailored specifically to suit this 
application A comnuter acauires full Raman spectra from the Raman 
spectrometer system, wherein data obtained from the spectra is processed 
using a multiple-dimension "least squares" routine 
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QP .« nE*rtm P TIQN OF THE DRAWINGS 

For a better understanding of the present invention and to show more 

clearly how it may be carried into effect, reference will now be made, by way 

ot example, to the accompanying drawings; 
5 Figure la illustrates a schematic representation of a Roman 

spectrometer system; 

Figure 1b illustrates a schematic representation cf a mirror 

arangement used m collaboration witn the cell illustrated in Figure la- 
Figure 2 illustrates a flow chart representation of a product 
10 diffeieiiLiahon algorithm used in determining product interfaces; 

Figure 3 illustrates a measured Raman spectrum of pure toluene 

obtained by the Raman spectrometer system Illustrated in Figure I; 

Figure 4 illustrates a measureo Raman spectrum f cr three grades of 

gasoline; 

15 Figure 5 Illustrates a graph of measured optimized correlation 

coefficients for five incoming pure compounds over a finite time period during 
the day, 

Figure 6 illustrates a graph of measured output voltage corresponding 
to transition interfaces between three grades of gasoline over a finite time 
20 period during the day, 

Figure 7 illustrates a graph of measured output voltage corresponding 
to transition interfaces between low sulphur, diesel oil, and jet fuel over a finite 
time period during the day; and 

Figure 8 illustrates a graph of measurea output vcliage corresponding 
25 to transition interfaces between low sulphur, furnace oil, and jet fuel over a 
finite time period uuiing the day. 

DETAILED DESCRIPTION OF THE INVENTION 

Figure la illustrates a Raman spectrometer system 10 for detecting 
30 interfaces between various petroteum products flowing through a pipeline or a 
suitable conduit structure It will also be appreciated that in accordance with 
the present invention, the product or material interfaces between other types 
of gas or liquid in a flow path can be detected. The spectrometer system 10 

wrnprieoo a high puwm wAuiLulior- »>wu;^ 12. whivh i..gluU*a «u high power <up 




to 500 mW), broad band diode laser device 14 with a centra! emission 
wavelength of 670 nm and an externa! cavity grating device 16 The diode 
laser 14 is locked intc single mode operation by means of the external cavity 
10, wherein the external cavity 16 provides frequency-selective optica! 
5 feedback. The diode laser 14 and externa! cavity 16 are operated to form a 
mode-locked external cavity laser that emits monochromatic red light with a 
spectra! resolution (FWHM linewidth) of 0.2 nm. The diode laser 14 also 
includes a thermoelectric cooling element (not shown in Figure 1a), a heat 
sink and fan, which provide* a means for maintaining wavelength stability and 

10 device integrity (avoiding device destruction). The Single mode output 
intensity is greater than or equal to the output from the laser whilst in 
muitimude operation when the laser is running at up to 07% capacity (as 
measured by the drive current). At drive currents above 67% of capacity, the 
single mode output intensity Is less than the Intensity measured during 

15 multimode operation. !n accordance with an alternative aspect of the present 
invention the excitation source may incorporate any other laser device 
capable of yeneiatiny wavelengths, which produce Raman spectra from the 
sample or products flowing in a carrier structure such as a pipeline. 
Depending on the material under analysis, rt may be necessary to use 

20 different excitation wavelengths for producing Raman spectra. Consequently, 
different diode laser sources may be utilized. 

The optical light signal (red light) is emitted from the excitation source 
12 and received by a lens device 18. The lens device 18 receives and focuses 
the light signal into the center of a high pressure, stainless steel sample eel! 

20 20. The cell 20 is equipped with four 1/2 inch diameter 4.75 mm thick AR 
coated fused sil;ca windows which allow the light signal (red light) to pass into 
and out of the cell 20. As illustrated in Figure 1b, a pair of yuld minus 202 
204, one below and one on the side of the coll respectively, are used to 
enhance the Raman signal intans'ty by effectively doubling the number of 

30 transitions (passes) the light signal makes through the sample. As illustrated 
in Figure 1b, a product under analysis, indicated by 206, flows into the ceil 20 
by means of a first opening 208. and flows out of the cell 20 by means of a 
second opening 210. As the product flows through the cell between the first 
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and second openings 208 210, the light signal, as indicated by 212, from the 
excitation source illuminates Ihe product, by making a first transition through 
the cell 20. The light signal, indicated by 212 (also referred to as an excitation 
wavelength) generates a Raman altered signal, wherein a first portion of 

5 the Raman signal is output from the cell along path 214 to the spectrometer 
device 26 (see Figure 1a) A second portion of the generated Raman signal is 
received by the first mi;rur 202 and reflected back along path 214 to the 
spectrometer device 26 (see Figure 1a) Consequently, the first mirror re- 
directs any generated Raman signal propagating away from the spectrometer 

10 device 26 back towards the spectrometer device 26, which increases the 
received Raman signal intensity The second mirror 204 reflects the light 
signal, indicated by 212, which passes through the flowing piuduut in the ceil 
back through the flowing product. Hence, by doubling the number of 
transitions the light signal, indicated by 212, makes through the cell, the 

15 Intensity o? the generated Raman scatteied signal s increased. Figure 1a 
illustrates the generated Raman scattering signal output from the ceil, as 
inoicated hy 22. The sample (product) flows continuously through the cell 20, 
which was pressure tested to 2000 psi. 



20 emitted from tne cell 20 and focussed onto an entrance slit 24 of a 
spectrometer 26 by means of a lens/filter assembly 2 8. The iens'filter 
assembly 28 comprises an optical filter 34, an input coupling lens 30 and an 
output coupling lens 32. The lens/filter assembly 26, removes any tesidual 
optical signals from the Raman signal (90% optical transmission above 680 

25 nm), which includes the 670nm light signal (red light) emitted from the 
excitation suuioe 12. The Raman signal 22 is coupled by the lens/filter 
assembly 28 to the entrance slit 24, where it is diffracted by a concave 
diffraction grating 36 within the spectrometer 26. The concave diffraction 
grating 36 directs the Raman signal onto the spectrometer s exit slit 38, 

30 wherein the exit slit 38 couples the Raman signal to an optical detection 
device 40 such as a photomultlpller tuce (PMT). The Raman siynal is 
detected by a photomultiplier tube 40 and sent to a computer device 42 for 
processing. A limited rotation electromag nehc drive (not shown in Figure 1a) 

icpctikively tu na the concave gristing 3€ through u jmail anglo under preciee 



As illustrated in Figure 1 



a, the Raman signal, indicated hy 22, is 




servo control. Turning Ifie yialiny 36 (one rotation or single scan) under servo 
control provides a wavelength scanning range cf 25nm, from 600-71 5nm. 
This corresponds; to Raman shifts ot 700-1220 cm" 1 . 

The Raman spectrometer system 10 is installed in-situ on a wording 
pipeline In a preferred embodiment, the excitation source 12, sample coll 20, 
spectrometer 26, and PM I 40 are housed in a NEMA rated explosion proof 
box with inside dimensions of 30 cm x 30 cm x 15 cm. This box is mounted 
outdoors: adjacent to the pipeline. The product is drawn through stainless 
steel tubing from the main pipeline, through the sample cell, and relumed to 
the main pipeline, providing a continuous flow through the cell A separate 
control unit 50 includes both the computer 42 and a power supply 46. and is 
located inside a nearby monitoring station. The power supply 46 provides 
power to all the main components of the Raman spectrometer system 10 (e.g 
computer, laser current drivar, cooling tan. etc ) A conductor meaium 44 such 
as a 20m length of single shielded co-axial cable transmits the Raman signal 
from the optical detection device 40 (photomultiplier tube) to the computer 42 
inside Thfi control unit 50. Once the computer 42 processes the received 
Raman signal, an indication signal consisting of a single voltage value scaled 
between 0 and AV is sent from the hack of the control unit 50 tc the end user. 
A change In the output voltage indicates a change in Lhe product flow'ng 
through the cell. More detailed data is aiso logged and stored on the computer 
42, and can be accessed for diagnostic purposes by including a computer 
monitor and keyboard within the control unit 50. 

The raw data based on the received Raman signal 22 is sent to the 
computer 42. This data consists of a continuous series of scans (2bnm 
range) each containing information capable or generating a Raman spectrum 
of optical intensity (response of the PMT optical detection device) as a 
function of wavelength (angle ot the diffraction grating) In accordance with 
the present invention, the data is averaged for approximately 90s t100 scans). 

The resulting incoming spectrum contains features which are 
representative of the composition of the product which is flowing through the 
cell. In accOT'ance with the present invention, a carefully constructed product 
differentiation algorithm is used to assign a signal value (an output voltage) 

UoacJ ui. Ll<* opKvoronoc of th<s nooming op«otrurr, ond/or tho extent to 
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which the incoming spectrum resembles previous incoming spectra, 
Consequently , the algorithm analyses the data corresponding to the averaged 
Raman scans, and determines the interface between different products 
flawing through the cell. In accordance with the present invention, the 
0 interface between petroleum products are determined, whereby the algorithm 
relies upon certain properties of refined petroleum products, which hava been 
established empirically in the laboratory during the development of the Raman 
spectrometer system 10. These properties are fairly baeic. and are expected 
to be widely applicable for most refined petroleum products world-wide. It will 

10 also be appreciated that this algorithm may be extended to analysing any 
substance that provides a good Raman signal. 

The Raman spectra of all grades of gasoline are dominated at Raman 
shifts between 700 and 1200 cm* 1 by features which can be attributed to a 
small group of compounds such as Toluene, nriho-, meia-, and psra-xylenes, 

15 Iso-octane. cyclohexane, and alkyl-cyclohexanes. Different grades of gasoline 
contain varying proportions of these compounds. For example, the toluene 
feature produces a Raman shift of 1001 cm 1 , and is always the most 
prominent. There is no universal relationship, which reliably correlates simple 
ratios of any of thase species to the octane rating of a given product. 

20 Any Raman features in spectra obtained from distillates such as diesel, 

low sulphur diesel, or furnace oil are obscured by a huge fluorescent signal. 
The relative strengths of the fluorescence signals from the different distillates 
were distinct and reproducible. Jet fuel fluoresced only weakly, but still 
showed little in the way of Raman features. Linear alkanes, which have very 

25 weak, broad Raman features between 700 and 120D cm" 1 , are present in 
greater abundance in distillates, including jet fuel. Aromatics and cyclic or 
branched alkanes, which have etrong Raman features, are undesirable in 
diesel fuels Occasional gasoline samples fluoresced weakly as they flowed 
through the pipeline, but not to the extent o f completely obscuring the Raman 

30 features. 

On the basis of these observations, the product differentiation 
algorithm was constructed along two paths. A flowchart mapping the decision 
making process is shown in Figure 2. 



Along the first path, strongly fluorescent products were assigned an 
output voltage in the range of 2.5-4 OV For the purpose of this disclosure, the 
exact value Is calculated based on the intensity of the fluorescence of the 
current sample, according to Equation [1]. For example, week!/ fluorescent jet 
5 fuel will have a very small value for F, and will produce an output value close 
to 2 6. Furnace oil, wnich can fluoresce so strongly the signal goes off scale, 
will have a value close to 4.0. The value of F in Equation [1] is constrained to 
a maximum ot H mcw = 2.0 Intermediate fluorescence levels associated with 
diesel cr low sulphur diesei, give intermediate (but distinct) output values. 

10 

V=2.5 + FM.4/F ma v [1] 
v value output to customer 

F ratio of fluorescence levels of incoming to stored reference spectre 
F ne > maximum fluorescence ratio before signal goes off scale 

15 

It will bo appreciated that a more detailed analysis of the incoming data 
is required to differentiate hptween gasoline products f-or this reason the 
algorithm must first determine whether the incoming product is gasoline or 
other distillates. Referring to Figure 2, in a step 56, the computer receives raw 

20 data corresponding to an acquired Raman spectrum (averaged spectrum 
based on the 90s scan time), Following this data acquisition process, in a step 
58.. the background fluorescence is subtracted from the acquired Raman 
spectrum data. In a step 60, a series of multi-dimensional least square 
calculations are carried out in order to fit the received incoming spectrum data 

25 to a set of data corresponding to the reference spectra already stored in the 
computer 42 (Figure la). The data corresponding to the reterence spectra are 
used in determining whether gasoline is being analyzed or other distillate 
products such as diesel, low sulphur diesel, jet fuel or furnace oil. The stored 
spectra (six) data are reference spectra of pure compounds, which were 

30 found to be important components of gasoline products. The stored reference 
spectra (data) are pure toluene, mefa-xylene, iso-octane, cyciohexane. and 
methylcyclohexane (five pure compounds), which were acquired in the 
laboratory and stored on the computer 42 (PC-104). Also, a sixth reference 
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spectrum (data) was added, which includes a characteristic curve from a 
strongly fluorescent compound In the step 60, data for each incoming 
spectrum is fit to the data corresponding to the mentioned reference spectra 
using a multi-dimensional least squares fit rcutine. The incoming product 

5 spectrum (data) is compared channel by channel (x values. Raman shift^m" 1 ) 
value) with a hypothetical "test" spectrum (data) which is a composite oT the 
six mentioned reference spectra (data). This tost spectrum data is built by 
applying a multiplier to ear:h reference spectrum, then taking the sum I He six 
multipliers are adjusted to minimize the difference between the incoming 

10 product spectrum (data) and this "test" or composite spectrum (data). The set 
of multipliers which minimizes the difference makes up a set of six optimized 
fit coefficients, one for each reference spectrum. The routine also calculates 
an individual correlation coefficient for each individual reference spectrum 
(data points) with respect to tne incoming product spectrum (data points). In 

15 the case where the incoming product contains substantial proportions of all or 
many of the individual compounds used to create the reference set, then none 
of these individual correlation coefficients will be especially close to one (i.e. 
the shape of the incoming product spectrum will not strongly resemble the 
shape of any individual reference spectrum). Hence for each incoming 

20 spectrum, this rcutine generates a set of six fit coefficients (the multipliers) 
and six individual correlation coefficients (a measure of the degree to whicn 
the incoming spectrum resembles a given reference spectrum), one pair for 
each stored reference spectrum. 

The routine also calculates a total correlation coefficient for the 

25 incuminy product spectrum (data) with respect to the optimized "test" or 
composite spectrum (data). If the product consists entirely or almost entirely 
of compounds which are included in the six reference spectra, the multi 
dimensional correlation coefficient will be very close to cne. This Lola! 
correlation coefficient gives a measure of how well the ircoming spectrum 

30 data matches tne composite reference spectrum data and used for diagnostic 
monitoring purposes. 

Based on the values of these individual correlation coefficients in a 
step 62, the algorithm tests fur the presence of strong Raman features in the 

incoming data If no Raman features ar© found, t^e output value ie calculated 



based on the fluorescence. In this case, the algorithm moves to a step 68, 
wh©roin the product ic determined to bo dioeel, jot fuel or furnace oil. If ctrong 

Raman features are present, the algorithm moves to a step 76, wherein the 
product is assumed to be gasoline. 
5 Subsequent steps 78, 80, 84, and 86 are for processing and detecting 

interfaces between varying grades of gasoline. In a step 78, data 
corresponding to the incoming spectrum is compared with data corresponding 

to 4 SfM-M:lriiMi which is i epieserilalive u! the piuduul lluwiny a lew minuleb 

previously. This previous product spectrum (data) is stored on the PC-104 
10 computer as a temporary or "moving" reference spectrum (data) It is updated 
every time the incoming spectrum is determined to contain strong Raman 
features (i.e. for every incoming spectrum, which indicates gasoline). In the 
step 78, the comparison between the incoming spectrum data and the moving 
reference spectrum data is performed using a single dimension least scuares 
15 fit The algorithm looks at the data corresponding to each spectrum (incoming 
spectrum and moving reference spectrum) point by point and compares the 
two. The difference between the corresponding signal intensities (the y 
values) for each channel number (the V values) can be quite large 
depending on how alike or unalike, with respect to both shape and signal 
20 intensity, the two spectra are. A constant multiplier is then applied to the 
stored "moving* reference spectrum data (i.e. multiplying all the "y" values 
corresponding to the stored mcving spectrum data by a given value). The 

multiplier that minimi coo the difference between the corresponding sicmd 
intensities (y values) is tne (best) fir coefficient. The correlation coefficient. 

25 normally denoted R PPl gives a measure of how good this fit is (i.e. how similar 
in shape the two spectra are). A perfect fit has a correlation coefficient of 
exactly one. Consequently, a fit coefficient and correlation coefficients are 
generated, which characterize the degree of similarity between the current, 
incoming spectrum, and the stored moving reference spectrum (the previous 

30 product spectrum), If the correlation coefficient Is above a defined cut-off. the 
product is deemed to be the same (no interface). In this case, in a step 82, the 
output value is naoilatetl according to equation [2a]. For this case, The output 
value is calculated from a previously calculated "moving" constant Cm defined 



by euuaUon [31 and the (small) difference between the incoming and previous 
product spectra (correlation coefficient). 

V = Cm + 1.0 + 5*<1-Rp P ) (if C f/ is unchanged or has increased! 
5 [2a] 

V * C M + 1.0 - 5*(1-Rpp) (if C w lias decreased) 

[2b] 

V value output to customer 
10 C w moving constant 

H pp correlation coefficient of incoming to previous product spectrum 

The third tarn is constrained to a maximum value of 0 5 to prevent V from 
going out of lanye. 

15 if the correlation coefficient ie below trie def nod cut-off this is taken to 

indicate 3 clear difference between the current and previous product spectra 
(an interface). In this case, in a step 80, the moving constant is re-assigned 
baaed on the chemical composition of the current incoming product according 
to Fquation [3] 

20 

Cm = A*(Rio + R<n-xyiw«) - B*R*e-ce [3] 

C M moving constant 

A. B empirically derived constants 

R. 0 , correlation coefficient of incoming spectrum to stored toluene reference 
25 spectrum 

FVxy*na correlation coefficient of incoming spectrum to stored meta-xylene 
reference spectrum 

R*o-c8 correlation coefficient of incoming spectrum to stored iso-ortane 
reference speclium 

30 In other embodiments of the present invention, whereby a product 

other than petroleum Is being analyzed, the equations ([1], [2] and [3]) which 
are used to determine the prndnct interface will vary In accordance with the 
type of product or comocund under analysis. For example, in accordance with 
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the present invention, the correlation coefficients defined in equation [3] relate 
to the particular incoming product (i.e. gasoline) being analyzed If another 
product is flowing through trie sample eel!, the correlation coefficients, will be 
Dased on other specific compounds stored as reference spectra in the 
5 computer, in the present embodiment, the correlation coefficients are based 
on the compounds found in gasoline, such as iso-octane, meta-xylene, and 
toluene. 

The constants A and B are empirically derived, and are used to scale 
the moving constant to give values between -0.5 and 1.0 for a wide range of 

10 gasoline samples. Once the moving constant has been updated, the output 
value is calculated using Equation [2a] or [2b] based on whether the recently 
calculated value of Cm is higher or lower than the previous C v value. If it is 
higher, then C M is substituted within equation [2a], and if lower, then C M is 
substituted within equation [2b]. In a step 84, the stored spectrum Is updated. 

15 whereby the present incoming spectrum data becomes the new moving 
reference spectrum (data), in a step 86, based on the different grade of 
gasoline detected al the interface, an output vullaye between 0-2. 5V is 
provided to the end user or customer (i.e, the output of equation [3]). The 
steps illustrated in Hgure 2 and described in the previous paragraphs 

20 constantly analyze incoming products (lowing in the cell. 

If the product ie determined to be dieeei, jet fuel or furnace oil, in a step 
70 the fluorescence level is calculated using equation [1). Depending on the 
product, an output value between 2.6-4.0V is generated and provided to the 
end user or customer, as defined by step 72. As previously discussed, each 

25 product (not gasoline} will have a relatively different fluorescence strength. 
Therefore, any changes in fluorescence signal strength can be correlated to 
determining a product interface. 

The system described above produced good quality Raman data ovei 
a wavelength span of 25 nm (Raman shifts from 700-1200 cm" 1 ). This shift 

30 region covers an aromatic ring breathing mode which gives a strong feature at 
990-1010 cm" 1 from mono-, nwta- or tri-substituted benzene (Important 
canetituomte of gasoline). It eloo eovero a ring vibration mode from pore 
substituted benzene (e.g. para-xylene. 810 and 830 cm '), and various modes 



from branched and cyclic alkanes (e.g cyclohexane ring breathing mode at 
802 err*', iso-octane f-butyl symmetric stretch at 745 cm" ) 

Figure 3 shows a spectrum, as defined by 100, of pure toluene 
ootained from this system 10 (Figure 1). The resolution is determined to be 
5 approximately 7.2cm* 1 FWHM for the 1001 cm' 1 peak, indicated at 102. Also, 
the pvriratlnn wavelength generated by the excitation source 12 is 6S6nm 
The ring breathing mode of toluene has a particularly strong Raman signal, 
and has proven to be a particularly useful benchmark feature for analyzing 
gasoline spectra. 

10 Refined gasoline is a mixture of hydrccarbons. Raman features from 

various aromatic compounds and cyclic or branched alkanes are visible in the 
spectra obtained from actual gasoline samples. Features from these 
compounds, which sir* impo'tflnt in defining ma nntane ra+ing nf a given fuel, 
are labeled in a set of gasoline spectra shown in Figure 4 These specua 

15 were obtained from three different grades of gasoline collected during fie;d 
trials at a pipeline pumping station located in north Toronto. The strong 
toluene features can be observed in all three spectra, as indicated at 104, 106 
and 100. The most obvious difference between ths three grades cf gasoline is 
the strength of the iso-octane peak at 745 cm 1 . Indicated at 110, which was 

20 clearest in the spectra of the premium grades of gasoline, indicated at 1 14 for 
the premium 91 grade, and at 112 for the US92 gasoline The iso-octane 
peak, as indicated at 116, in the spectrum of regular unleaded gasoline 
(lowest octane rating) was even loss distinct. A closer examination reveals 
other differences in the three spectra, visible predominantly between 700 and 

25 775 cm" 1 . For example, the /nete-xylene peak at 725 cm" 1 , as indicated at 120, 
is most prominent in the spectrum of US92 premium grade gasoline, indicated 
at 122. while regular un'eadao fuel has stronger cyclohexane and alkyi- 
cyclohexane features. 

The chemical differences between regular vs the premium grades of 

30 gasoline as indicated by the gasoline samples illustrated in Figure 4, appear 
distinguishable. Aromatics and branched alkanes such as iso octane are 
recognized octane enhancers It is thus reasonable that these compounds 
should feature more prominently in the spectra of the premium grades. The 

epactra shown in Figure 4 clearly establishes that there are diffo-oncos 
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between gasoline samples, and that these differences can be detected 
spcctroscopically. 

Hor the products measured over the course of tvvc lengthy field trials, 
mere was a tendency for the premium grades tc have higher proportions of 
6 iso-octane relative to toluene than the lower (regular) grade. However, the 
ratios from some samples were Inverted (I.e. nigh levels of iso-octane in some 
regular unleaded samples, low levels in some premium US82 fuels) Also 
there are some instances in which no clear change in the iso-octane to 
aiumaliu (toluene) ratio can be observed, even when a transition is krown to 
10 have occurred. As demonstrated, these •missed" transitions are not due to a 
lack of sensitivity in the spectrometer Rather, they are not detected because 
many different compounds, not just iso-octane and toluene, affect the overall 
octane rating of a given fuel. 

In accordance with the present invention, the end user will detect the 
15 required transitions (product interfaces) by a change in a single value (e.g. 
voltage), which is output as a voltage and integrated within the user's pipeline 
rnuniUning station. The noise level in the output value has to be sufficiently 
low relative to the magnitude of the expected changes, eo that an interface 
could be reliably identified within one or two measurement periods. 
20 In order to overcome the mentioned issues with detecting the gasoline 

transitions, the multi-dimensional least squares fitting is used. The set of six fit 
coefficients and correlation coefficients, which are calculated basea on the 
corresponding six stored reference spectra, provide a quantitative basis for 
differentiating between gasoline samples. The magnitudes of the fit 
25 coefficients depend on the strength of the Raman signal, and thus are 
sensitive to external influences such as the laser intensity, or dirt particles or 
accumulation in the cell, which may vary over time Hnwever, the rwreiatif»n 
coefficients give a measure of the degree to which the incoming spectrum 
resembles a given reference compound, which proves to be a more reliable 
30 indicator of a pmritic.t r.hanrjfi. This is demonstrated in Figure 5, which shows 
Hie optimized correlation coefficients of the incoming product spectrum 
calculated with respect to five pure compounds. Three product interfaces 
occurred In a 24-hour period. As illustrated in Figure 5, no singie correlation 

coefficient responded reliably at ever/ product interface, f" jrthermcre, the 



magnitude of the change* at some interfaces were small relative to the 
degree of scatter giving an unsatisfactory signal to noi6e ratio. These 
problems are overcome by calculating a single output value by using 
equations [2a, 2b and 3] in accordance wiih the present invention. As 
5 previously described, the output value ie derived from the correlation 
coefficients ot the pure compound reference spectra, and the correlation 
coefficient of the current (incoming) spectrum to a stored spectrum of the 
product which had been flowing in the previous few rrinutes. Figure 6 shows 
this calculated output (In real or near real-time) value plotted as function 
10 gasoline flow over a period of 24-hours, The output value, defined by 140, 
reliably detects gasoline interfaces, as indicated at 142. 144, 146, and gives a 

yyuU fciyiieil lu mjiae iciliu. Tui ca ample, the mlir« r «oc, inJicALed by 14£, 

between regular unleaded and US 92 is defined by a relatively well defined 
voltage step, as indicated at 150. 

10 A further advantage of this method for detecting gasoline interfaces is 

that the algorithm was designed to exaggerate the difference in output values 
right at the interface. This is useful in cases wheie two dilTeierit grades of 
gasoline have chemical compositions, which are very similar with respect to 
the compounds to which the Raman spectrometer is sensitive A good 

20 example of this exaggerated change is seen In Figure 5, where a transition at 
the interface, ae indicated at 146, between premium 91 and regu'ar unieaded 
occurs at approximately 20:30 hours. At the transition, the output value dips 
sharply from a value of 1 .82 V to a minimum of 1.22 V, before rising again to 
level off at 1.56 V. The fnal values are related to the chemical composition of 

25 the product, while the size of the dip Is related to the degree of difference 
between the incoming (current) product spectrum and the stored (previous) 
product spectrum. Two chemically similar gasoline products will give very 
similar final output values. However, the transition is still detected due to the 
small (but measurable) change in the product spectrum, which occurs at the 

30 interface. 

As described in the flow chart of Fiyuie 2, a different strategy 
(fluorescence signal strength) is used to detect product interfaces between 
tne distillates This successfully differentiates between jet fuel, diesel oil, iow 

oulphur dicool, Orid fumaoe oil (not ohovvn) ao illuoiratod in iho firot fiold ir&i\ 



results shown in Figure 7. Figure 7 provides a plot uf output voltage, which is 
calculated using equation [1], against time of day. Figure 7 shows a quick set 
of transitions from low sulphur riiesel to diesei oil, as indicated at 150 diesei 
oil to jet fuel, as indicated at 152 followed by a return to diesei oil, indicated at 

5 154, then back to low sulphur diesei, as indicated at 156. No product was 
flowing through the pipeline between 1)8 3b and 11:15 hours. Particularly 
noteworthy among this set of transitions are the clear differences between 
diesei oil and lew sulphur diesei, as these products can be difficult to 
distinguish by previously existing technologies Transitions from a gasoline to 

10 a distillate or vice versa (not shown) are easily detected by this, and other, 
techniques 

Data from a second field trial is shown in Figure 8, wherein fie 
interface transitions between jet fuel, low sulphur diesei oil, and furnace oil are 
shown The transitions between iow sulphur diesei oil ana furnace oil, as 

15 indicated at 160, between furnace oil and jet fuel, indicated at 162, jet fuel and 
furnace oil, indicated at 164, and furnace oil and low sulphur diesei oil, 
indicated by 166, are clearly shown. Consequently, In accordance with the 
present invention, the Raman spectrometer system provides a clear indication 
of interfaces between petroleum products flowing in a pipeline The system 

20 detects the Interfaces between both various distillates and various grades of 
gasoline products by providing the end user or customer with a single 
indication value. 

Obtaining good quality Raman dala ib a Key inyiedienl foi successfully 
detecting interfaces in a flowing product stream. The Raman spectra acquired 

25 by the Rgman spertrnmerer system has excellent resolution and a gone 
signal to noise ratio. The resolution of a Raman spectrum in general may be 
limited by the resolution (degree of monochromatically) of the excitation 
source, or by elements within the spectrometer. By using an external cavity to 
mode lock the diode iaser, an excellent degree of monochromatically is 

30 achieved. If the excitation source is the limiting factor defining the resolution o* 
the instrument, the observed excitation source linewidth (FWHM) of 0.2 nm 
will translate into a Raman shift of 4.5 cm*". The observed resolution of 7.2 
cm* 1 (at the hWHM) is actually limited by the nnmhinerJ diffraction 



yrating/fccuswng element inside the spectrometer, which has a dispersion of 
1 mm/5 nm. 

The nigh output power of The laser, high transmission tnrough the 
single element spectrometer (low transmission losses) and low dark current 

5 levet of the PMT detector all contributed to establishing a good signal to noise 
ratio. This means the scan time can be decreased to approximately 90s. This 
sampling interval was 3hort enough to detect product interfaces in real or 
near-real time. Ideally, the interface between products is sharply defined In 
practice, souk* mixing occurs in the pipeline, depending on the level of activity 

10 (product flow rates} through the line. Real interfaces can thus span a period of 
up tn 9 minutes a time sinn If Inaptly rjrRRter than the sampling period 

The spectrometer as designed is compact and rugged, idea! for in-situ 
use in hostile environments. In this application, the spectrometer unit was 
located outdoors, adjacent to the gasoline pipeline. The sampling method is 

15 noninvasive and non-destructive (i.e. not a drop of product was consumed in 
the testing). This feature is also attractive for situations In which contamination 
may be an issue^ for example in the two extreme cases in which the sample i3 
either designed for human consumption or is highly toxic 

Comparing the incoming (current) product spectrum with a stored 

20 spectrum of product, which had been flowing a few minutes previously, 
proved tc he the most sensitive method for detecting gasoline to gasoline 

product into Tow*. Tins algorithm w** de*lgn*d Lw produce an eActyyciated 

difference in output values right at an interface, then stabilize over a period of 
a few minutes to an output value related to the bulk chemical composition of 

25 the current product. This feature is invaluable for distinguishing between 
gasoline of very similar compositions. 

The software in general was designed to be lubusl. Distinguishing 
between gasoline on the basic cf correlation, rather than fit coeffjeiente 
RliminnTAB noise" due to changes in the Kaman signal intensity. The software 

30 was also designed to compensate for small gradual or sudden changes in the 
position (channel number) of the major Raman. The software corrects for this 
by "tracking" the toluene ring breathing mode, which has a Raman shift of 
1001 cm" 1 . In the region scanned, this peak always appears as the strongest 
feature in tho Raman epootra of goeolm©. The software routine find* this peak 



in the incoming product spectrum and compares its position with that of the 
analogous feature in the stored reference spectrum of toluene If the position 
(1001 cm 1 ) of the Toluene peak (incoming spectrum) shifts, the shift is 

uuncUcu ao that the maximum of the telucne peak in tho incoming opootrum 
6 occurs at the same channel number (1001 cm' 1 ) as in the reference spectrum 
The fining routines are then performed as normal on the uui reeled (shifted) 
incoming data. Such shifts could arise from a ch.ango in the behavior o f the 
scanner motor, a sudden shift (mode hop) in the excitation wavelength, or a 
yiadual shift in the excitation wavelength due to temperature induced changes 

1 0 in the alignment of th^ external cavity. 

In accordance with the present invention, the Raman spectrometer 
system described herein, is not limited to specifically detecting interfaces 
between petroleum products. The system and product differentiation algorithm 
can be utilized to detect product interfaces in a variety of yas oi liquid 

15 products, by means of accordingly varying the excitation source wavelength 
and other wavelength dependent components (e.g. gratings) In the system 



